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I Abstract. We derive the fluctuating hydrodynamic equation for the number and 

momentum densities exactly from the underdamped Langevin equation. This 
derivation is an extension of the Kawasaki-Dean formula in underdamped case. The 
steady state probability distribution of the number and momentum densities field can 
'"^ . be expressed by the kinetic and potential energies. In the masslcss limit, the obtained 

C I fluctuating hydrodynamic equation reduces to the Kawasaki-Dean equation. Moreover, 

the derived equation corresponds to the field equation derived from the canonical 
equation when the friction coefficient is zero. 
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1. Introduction 



o 

OO . Field equation is widely employed in the studies on colloidal or liquid dynamics. In 
^ ' the study on colloidal dynamics, some researchers have applied averaged density field 
dynamics, which is called the time-dependent density functional method [1-7]. This 
^ I method has also been successfully employed to study various phenomena observed in 
I the field of liquid dynamics, such as solvation [7-13], transport phenomena [14], and slow 
relaxation in supercooled liquids [15,16]. Besides the average density field dynamics, 
other researchers have also developed theoretical expressions describing momentum 
density fields [17, 18]. 

As compared to the direct calculation of particle dynamics, field description is more 
useful for theoretical studies. This is because we can estimate many physical parameters, 
such as transport coefficients from the correlation functions of field variables. Thus, 
by using the field description, many researchers have formulated approximations for 
the estimation of the physical parameters. For example, the mode-coupling theory, 
which is known as a useful tool for approximation of the transport coefficients, has been 
formulated by using field description [19-22]. 
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While the field description is useful for theoretical calculations, its correspondence 
with the particle description is not clear. Therefore, the derivation of the field description 
from the particle description is the fundamental problem in the studies on colloidal and 
liquid dynamics. When field variables are not averaged. Dean has derived the field 
equation from the overdamped Langevin model [23]. In a colloidal system, the time- 
dependent density functional method can be applied for the derivation of field equations 
from the overdamped Langevin equation by averaging the density field on the basis of 
some assumptions [1-3]. Recently, a method has been developed to derive field equations 
from the Liouville equation describing liquid dynamics by using the projection operator 
method [24,25]. In most cases, the derivation of the field description from the particle 
description requires some approximations. 

Very few studied have been carried out on the derivation of the field description 
in the nonlinear and underdamped cases. In these cases, the inertial effect has to 
be considered. In liquid dynamics, linear generalized Langevin equations including 
momentum density have been derived for the field variables of a homogeneous system 
[19]. Linear generalized linear Langevin equations for an inhomogeneous system have 
also been developed [17, 18]. However, nonlinear equations in the filed description have 
not been derived. Therefore, in underdamped cases, phenomenological models have 
often been employed [26]. 

In Ref. [23], the evolution equation of the density field is derived from the 
overdamped Langevin equation representing the particles interacting via the pairwise 
potential. The derived equation is called as the 'Kawasaki-Dean formula'. In Ref. [23], 
a closed evolution equation for the density field is exactly derived by using Ito's formula 
[27], while the evolution equation is approximately derived using other field models. The 
steady-state probability distribution of the density field for the overdamped Langevin 
model is represented by the bare pairwise potential term and the entropy term. In 
contrast to the overdamped case, there are no exact derivations of the evolution equation 
for the field variables in the underdamped cases. Here, a question arises whether we can 
extend the Kawasaki-Dean formula to the underdamped Langevin equation. The exact 
derivation of the closed evolution equation for field variables is main issue of this paper. 

In Sec. 2, we derive the closed evolution equation for the number density and 
the momentum density field using the underdamped Langevin model. For a system 
without dissipation, the closed evolution equation corresponds to the field equation for 
a Hamiltonian system. In Sec. 3, we discuss the properties of the derived evolution 
equation. In the Sec. 3.1, we calculate the steady state probability distribution 
functional of the evolution equation by using a functional Fokker-Planck equation. In 
Sec. 3.2, we derive the Kawasaki-Dean equation from the evolution equation derived in 
Sec. 2 to check the consistency between our model and other models. Sec. 4 presents 
the concluding remarks. 
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2. Derivation of the nonlinear fluctuating hydrodynamic equation from 
underdamped Langevin equation 

We study Brownian particles suspended in a three-dimensional solvent at temperature 
T. The motion of the i-th Brownian particle is represented by its position Xi and 
momentum p^, where i = 1,2, N and Xi G [0, L] x [0, L] x [0, L]. We express the a-th 
component of Xi as x", where a = 1,2, and 3. That is, Xi = {xl,xf,x^). The Brownian 
particles interact via the pairwise potential V{x). Each Brownian particle has the same 
mass m. The motion of the i-th Brownian particle is described by the underdamped 
Langevin equation as 

dx^^Pi 
dt m' 

^-^-^■^^-^^^ 

where U{{xi}fL^) denotes the total potential energy defined as 



i=l j=lj^i 

The coefficient 7 is the friction constant and Ri{t) is the zero- mean Gaussian white 
noise satisfying 

Rt{t)R^{t')'^ = 25,,5^pSit-t'), (4) 

where (■) represents the average value of Ri{t). 

First, as described in Ref. [23], we introduce the density field p{x,t) given as 

N 

p{x,t) = Y,^{x-x,{t)). (5) 

i=l 

To obtain the closed evolution equation of the fields in the underdamped system, we 
also introduce the momentum density fields g{x,t) defined as 

N 

g{x,t)=Y,P^{m^-X^{t))■ (6) 

i=l 

One can exactly derive the closed evolution equation of the number density field and 
momentum density fields defined by equations (5) and (6). 

Using these definitions, we derive the evolution equation using Ito's formula. 
Expanding the stochastic variable p{x, t) defined in equation (5), we obtain the evolution 
equation for the number density as 

dp{x,t) ( 9{x,t) \ 

dt \ m )' ^ ' 

where we have used equations (1) and (6). This equation represents the continuous 
equation for the density field. Similarly, using equations (2) and (6), we obtain the 
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evolution equation for the momentum density 

ot m 

- p{x, t) j dx^^-^piy, t) ^yi^ (8) 

where we have used Einstein's summation convention whenever a subscript is repeated 
in a term. Here, ^'^{x,t) and M"^ are defined as 

N 

C{x,t) ^ Yl VlTRtrnx - X,), (9) 

i=l 

M"''{x,t)= yPl^Mi^S{x-x,{t)). (10) 

i=l 

Trivially, the average of $, is zero from equation (9). Further, ^ is a multiplicative 
noise: time correlation depends on the instantaneous density fields. The noise term in 
equation (9) is rewritten in the form 

C{x,t) = ./T-P{x,t)TC^{x,t), (11) 
where T°'^{x,t) is defined by 

T''^{x,t)=^p{x,t)6afs (12) 
and C is the space-time Gaussian white noise satisfying 

(C(x, t)Cf{x', t')) = 26o.p6{x - x')6{t - t'). (13) 

To obtain the closed evolution equation, we have to make the following assumption 
for the trajectory of the positions of a particle {Xi{t)}f^^: 

S{xi{t) - xj{t)) = 6{xi{t) - xj{t))6i,, (14) 

where Sij is the Kronecker delta. Note that our aim is to construct a map from 
the trajectory of the position and momentum of the particles {xi{t) , p^{t)}'^Li to the 
trajectory of density and momentum density fields [pt,gt] = {p{x,t), g{x,t)}x- Then, 
equation (14) is satisfied when no two particles occupy the same position simultaneously 
in the mapping. Such an assumption is valid if particles interact through via a repulsive 
pairwise potential and a discretization of space, which is discussed in Appendix A. 
Equation (10) is formally rewritten in the form 



m 

i=l 

N N 



mp[x,i) .^^ 

The infinite form included in (15) is explained in Appendix A. In the second step, we 
have used the definition of the density field given by the equation (5). Then, by using 
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equation (14), M°^(a;,t) is represented only by [pt,gf] in the form 

^ N N 

^ g-ix,t)g''ix,t) ^^g^ 
mp{x,t) ' 

where we have used the definition of the momentum density given by equation (6) in 
the last step. 

Finally, substituting equation (16) into equation (8), we obtain the evolution 
equation for the momentum density as follows: 

- o(x t)^ f^lMpy] _ ^ ( 9''MAx,t) \ 
^^''dx^\5p{x,t)) dxP\ mp{x,t) J' 

where we use the abbreviation for the functional derivative as 
6Hv[p] 6Hv[v] 



5p{x,t) 6(p{x) 
6HK[p,g] 6Hk[^,iI)] 



ip(x)=p{x,t) 



6g°'{x,t) 6lp°'{x) ^(x)=p{x,t),ip{x)=g{x,t) 

This abbreviation is used hereinafter. Here, -ffyfv?] and Hxlf, V'] ^ire functionals for the 
functions ip{x) and if^{x), respectively, and are defined as 

HvM =^ J dx j dyV{x - y){Lp{x)Lp{y) - 5{x - y)ip{x)), (19) 

Clearly, the functionals in equations (19) and (20) correspond to the internal energy and 
the kinetic energy of the system, respectively. Equations (7) and (17) are the desired 
nonlinear fiuctuating hydrodynamic equations. 

The first term on the right-hand side of equation (17) leads to the decay of 
momentum. The dissipative matrix F"^ in equation (17) depends on p{x). This 
feature, which is a characteristic of the Brownian particle system, is in contrast to 
features of the Navier-Stokes equation. The dissipative matrix in the Navier-Stokes 
equation is given by the combination of the gradient and the shear and bulk viscosities 
[21,26]. The fiuctuation-dissipation relation of the second kind is satisfied by the first 
and second terms on the right-hand side of equation (17). That is, the dissipative matrix 
is consistent with the noise coefficient. 

The third and fourth terms on the right-hand side of equation (17) represent the 
conservative fiows. The fiow represented by the fourth term is caused by the momentum 
transfer. Further, the fiow represented by the third term is caused by the gradient of 
the functional derivative of the Hamiltonian including the bare potential V{x — y) (or 
Hv[p]) in equation (19). This is in contrast to many field models including the chemical 
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potential or free energy. The bare potential is obtained by the exact derivation from 
the overdamped Langevin model [23] . Thus, the present result shows that the inclusion 
of the bare potential is general consequence of the exact derivation without any coarse 
graining. 

The Hamiltonian in equation (17) does not include the entropy terms, which 
are included in the overdamped evolution equation for the fields [23] or in the 
phenomenological model in the underdamped case [7, 28] . One can drive the entropy 
terms for the Brownian particle system from the momentum transfer term when the 
overdamped limit is considered in equation (17) (Sec. 3.2). The entropy terms for liquid 
dynamics also originate from the momentum transfer term in the Liouville equation [25]. 
Those results indicates that the entropy term in the evolution equations is eliminated by 
explicitly treating the momentum transfer term from the point of view of the derivation 
from a microscopic model. 

From equations (7) and (17), we also obtain the closed evolution equation of the 
density and momentum density from the canonical equation. Equations (1) and (2) 
reduced to the canonical equation when 7 = 0. Therefore, by substituting 7 = into 
equations (7) and (17), we obtain 



dp{x, 
dt 


t) _ 


- V 


'9{x 

m 


t)] 






dg{x 
dt 


t) _ 


- pix, 


t)V 


iSHvip]] 
_Sp{x,t)_ 


- V ■ 



'g{x,t)g{x,t) 



(21) 
(22) 



mp{x, t) 

The evolution equations (21) and (22) contain the following five conserved quantities: 
the total energy, the total number, and total momentums. These quantities are defined 

as 



H[p,g]^Hv[p]+HK[p.g] 
N[p] = [ dxp{x), 



P[g] = J dxg{x). 



(23) 
(24) 

(25) 



The conservation law for the total energy functional is derived as follows: 



dH[pt,gt 
dt 



dx 



dp{x,t) 
dt 



5Hv[p] g\x,t) 



6p{x,t) 2mp'^{x,t) 



dg^{x,t) g^{x,t) 



+ / dx^^^^^p^^^^p^. (26) 
J dt mp{x,t) ^ ' 

By substituting equations (21) and (22) into equation (26) and integrating by parts 

several times, we have 

dH[pt,g^] 



dt 



dxV 



g{x,t) f 5Hv[p] 6HK[p,g] 



m 



(27) 



6p{x,t) 5p{x,t) 

The left-hand side of equation (27) is equal to zero from the divergence theorem. In 
addition, the conservation law for the total number can be easily checked from equation 
(21). 
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The conservation law for the total momentum is also proved as follows: 



dt 



dxp{x, t)V 



SHv[p] 



dxV ■ 



g{x,t)g{x,t) 
mp{x, t) 



6p{x,t)_ 

The first term on the right-hand side of equation (28) vanishes by the action-reaction law. 
The second term on the right-hand side of equation (28) vanishes from the divergence 
theorem. Therefore, the total momentums are conserved. Note that we have obtained 
the conservation law directly from the continuous model given by equations (21) and 
(22) without using the canonical equations (1) and (2). 

Equations (21) and (22) are similar to the Euler equation in fluid mechanics [29]. In 
these equations, the number and the momentum are conserved, and the advection term 
is present. However, there are some differences between them, which will be discussed 
in Sec. 4. 



3. Properties of the nonlinear fluctuating hydrodynamic equation 

In this section, we discuss some aspects of the closed stochastic evolution equations (7) 
and (17) along with the Hamiltonians (19) and (20) and the noise given by equation 
(13). 

3.1. Derivation of Fokker-Planck equation for underdamped fluctuating hydrodynamic 
equation 

In this subsection, we calculate the steady-state probability distribution functional 
for the number and momentum density fields from the derived stochastic evolution 
equations (7) and (17). We first derive the Fokker-Planck equation for these field 
variables by using a standard procedure. Then, we obtain the steady-state probability 
distribution functional as a stationary solution for the Fokker-Planck equation. In 
this subsection, a time-dependent function f{x,t) is denoted by ft{x) using standard 
notations for a stochastic process. 

The probability density distribution functional is defined as 

Pi[p,glt) = {S[p-pt]5[g-g,]), (29) 

where (■) represents the average over {(^^{x)}x. is a delta functional defined as 

5[p - pt]6[g -g,]^ll6 (p(x) - p,{x)) 5 {g{x) - g,{x)) . (30) 

The evolution equation for the density field given by (7) is rewritten in the form 

dpt{x) = —V ■ g^{x)dt/m. (31) 

In addition, the evolution equations for the momentum density fields given by (17) is 
rewritten in the form 

dg,{x) = G{pt{x),g,{x))dt+ ^/jpt{x)Tdri,{x). (32) 
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Here, G is defined as 



and rj'lix) satisfies 



5HK[p,g] , d 



5Hy[p\ 

5p{x) 



d 



g"{x)g^{x 
mp[x) 



d'nUx)d7]tix') = 25a(i5{x - x')dt. 

Here, drj"{x)d'r]^,{x') is equal to zero in the case olt ^t' . 

To obtain the Fokker-Planck equation, we apply Ito's formula in (30) as follows 

d{5[p- Pt]5[g - Qt]} 

5{5[p-pMg-g,]} 



(33) 
(34) 



+ 



dxdpt{x) 
dxdg^{x) ■ 
dxdg^{x 



5pt{x) 
6 {6[p - pt]6[g - g,]} 



6 



6gt{x) 



dx'dg^{x 



6 {6[p - pt]S[g - gt]} 



5g,{x') 



(35) 



Here, we have defined the functional derivative as 
6{6[p-p,]5[g-g,]] 



dxdpt{x) 
lim > dptixj) 



5pt{x) 
d 



dpt{xi] 



X\5{pt{xr) - p{xr))6{gt{xi') - g{xr)), 



(36) 



where / and /' are the indices of the discretized space coordinate with volume | Aa;| and 
xi is the discretized position. By substituting the evolution equations (31) and (32) 
into equation (35), we obtain 

d{^[p- pMq- Qt]} 



— / dx'V ■ gf{x)dt/m- 



6pt{x) 



dx {G{pt{x),g^{x))dt+ ^/-fTpt{x)d'nt{x] 
^ J dx^-iTpt{x)dr)t{x) 



8 {8[p - p,]5[g - g,]] 



6gt{x) 



6gt{x) 



dx' y jTpt {x')drj^ {x 



/^ ^ {^[p - Ptl^ia - 9t]} 
5g,{x') 



(37) 



Changing the index of the derivative of the delta functional from [pt.Qti to [p, g], we 
obtain 

d{^[p~ Pt]5[9 - at]} 

dx-^^ [V ■ g{x)dt/m {6[p - pt]6[g - g^]}] 



Sp{x) 



— / dx- 



d 



+ dx 



5g{x) 
6 



6g{x) 5g{x] 



G{p{x),g{x))dt + y/^Tpt{x)dr),{x)^ {6[p - pt]6[g - g,]} 
^ [^Tp{x)dt{d[p-p,]6[g-g,]}]. 



(38) 
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Here, we have integrated the third term on the right-hand side of equation (38) with 
respect to x' using equation (34). By substituting equation (33) into equation (38) and 
averaging (38), we obtain the Fokker-Planck equation for [p,g] as 

^^^I|^ = £([p,^])P([p,y],t). (39) 

Here, the operator £ is a hnear operator defined as 

S — f g{x 



^i[P,9]) = J dx 



6 

+ 



5g{x) 



Sp{x) \ m 



Sp{x) J \ mp(x 



6g{x) V ^9{x) 6g{x) J] 

where V^'^i^x) = 'yp{x)6a/3. We obtain the steady-state probabihty distribution 
functional Peq[p,g] as a stationary solution of equation (39) given by 

1 / Hv[p]+H K[p,g] 
Z'^^P [ T 

where Z is a normalization constant determined by 



^cq[p, g\ = - exp ( ) , (41) 



VpVgP,q[p,g] = l. (42) 

Note that Peq[p, g] does not include the entropy terms although the Hamiltonian in 
the steady-state distribution obtained using the overdamped field model includes them 
[23]. Further, the steady-state distribution functional Pcq[p,g] has same form as the 
classical fluid [7]. 

3. 2. Massless limit of the underdamped fluctuating hydrodynamics 

In this subsection, from equations (7) and (17), we derive the overdamped fluctuating 
hydrodynamic equation for Brownian particle systems in the massless limit. In the 
massless limit, the obtained equation is the so-called Kawasaki-Dean formula. Note that 
the Kawasaki-Dean formula is derived from the overdamped Langevin equation, which is 
obtained from the underdamped Langevin equation in the massless limit. Therefore, the 
derivation of Kawasaki-Dean formula from our equations leads to a consistency between 
our equations and these equations. Moreover, the derivation given in this section is 
useful for understanding similar studies carried out in the past [7, 25] . 
Using equations (5) and (14), we obtain the identity 

p{x, t)p{x, t) = 6{x — x)p{x, t). (43) 

Here, the right-hand side of the identity has an infinite value, which can be justified by 
the proper interpretation of discretization discussed in Appendix A. Equation (43) is 
satisfied only if the density is defined as the sum of delta functions. Therefore if the 
density is defined as a continuous function, equation (43) is not satisfied. 
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By using r = rn/'-f, the evolution equations of the density field and momentum 
density field are rewritten in the form 
dp{x,t) 1 



dt 

dg°'{x, t) 
dt 



-V-g(a;,t) 

g'^ix.t) 



(44) 



_d_ 



d 



g^{x,t)g^{x,t) 



6p{x,t)_ 

+ y/jTp{x,t)Cix,t). 



(45) 



T'jp{x,t) 

The parameter r is the relaxation time for the density field and is constant for a 
given system. We focus on the time evolution of the density field whose time resolution 
At is significantly larger than r. Then, we define 



p{^, tn) 



lim p(x,tr, 

r/At-^O 



(46) 



where tn = nAt. After taking the limit of r and evaluating the equation, we take the 
continuous limit At — and represent the time evolution of the coarse-grained density 
field p as follows: 



dp{x,t) 

at 



= lim 

At^O 



lim 

At^O 



p{x,tn + At) - p{x,tn) 
At 



lim 

T-»0 



p{x, t + At) — p{x, t) 
At 



(47) 



In the derivation of the coarse-grained evolution equation, we have used the 
asymptotic formula 



lim 



t' 



t' -{t'~t)/r 

dt- A{t) = A{t') for t' > t'_. 

r 



(48) 



That is because a term in the integrand is used in the definition of the delta function 



5(t - t') = hm 



-\t-t'\/r 



(49) 



-0 2r 

Note that t' is the upper limit of the integral in equation (48). 

By integrating (44) with respect to time, we obtain the difference p{x,t + At) — 
p{x,t) in equation (47) as follows: 

1 /•*+^* 1 

p(x,t + At) -p(a;,t) = V-/ dt'-g{x,t'). (50) 

1 Jt ^ 

Next, we consider a system with t ^ r. By using equation (45), g{x,t') is formally 
solved as 

-g{x, t') = -g{x, 0)e-*'/^ + r{x, t') + H(a;, t') + n(a;, t'). 



T T 

Here, we have introduced the quantities 



T-{x,t')^-ll^ dse-^''-^y^p{x,s)^ 



6p{x, s) 



(51) 



(52) 
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rs=t' 



E'^{x,t') ^ 



e-(*'-)/V7Tp(aj,s)ci<(a;), 



s=0 



d 



g'^{x,s)g^{x,s) 
T-fp{x, s) 



(53) 



(54) 



The terms H, T, and 11 correspond to noise, drift, and advection terms respectively. 
Substituting equation (51) into equation (50), we obtain 

-t+At 



-f{p{x,t + At)- p{x,t)) 



dt' (T(a;, t') + H(a;, t') + n(a;, t')) 



(55) 



Using equations (52), (53), (54), and (55), we evaluate the right-hand side of 
equation (47) as follows. First, we substitute equations (52), (53), and (54) recursively 
into the left-hand side of equation (55). Next, taking the limit r ^ 0, we evaluate 
it to the first order of At. Then, taking the continuous limit At 0, we obtain the 
right-hand side of equation (47). 

First, we integrate of T. By using equation (48), equation (52) is evaluated as 



limT"(a;,t') 



-p{x,t 



5p{x,t') 



Then, the integral of T{x,t') in the limit r — > is evaluated as 



lim 



t+At 



dt'T'^ixA') 



p{x,t) 



_d_ 



6p{x,t) 



At. 



(56) 



(57) 



Next, to integrate H, we calculate the correlation for these variables in the case of 
At ^ T. The product 'E°'{x,ti)'E^{x',t2) is calculated as 



'(x,ti)H^(x',t2) 



si=0 Js2=0 



--lT^p{x,s^)p{x\s2)dr^l ix)dr]^^ (x') 



2-fTS{x - x')6''''e 



'\^^l3f.-\ti-t2\/r 



min[ti,t2] g-2(min[ti,t2]-s)/T 

ds p{x,s). (58) 



Here, we have used Ito calculus (34) and the identity ti + 12 = |ti — t2| + 2min[ti,t2]. 
By integrating equation (58) with respect to time ti and t2, we represent the product 
of integrations of S as 



t+At 



dti 



t+At 



cit2H"(a;,ti)H^(a;',t2) = 2^T6^p6{x - x') 



t+At 



dti 



t+At 



dto 



-\tl-t2\/T 



X 



min[ti,t2] 



ds- 



-2(min[ti,t2]-s)/T 



-pix,s). (59) 



Taking the limit r — in equation (59), we obtain 



/t+At Pt+Llt 
dti cit2S"(a;, ti)H^(a;', t2) = 2-fT5a(i6{x - x')p{x, t)At + o(At). (60) 

Comparing equation (60) with equations (11), (12), and (13), we find that the time 
average of H(a;,t) coincides with that of ^{x,t) when p is replaced with p. 



ft+At 



lim 
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Next, we integrate of 11. By substituting equation (51) into equation (54) 
recursively, we integrate of Il{x, t') as follows: 



t+At 



dt'U'^ixA') 



d 



d 



t+At 



dt' / dse-^''-'^'^ 



7p(a3, s) 



T 



t+At 



dt' / dse 



-fp{x, s) 



+ o(At). (61) 



In the second step, we have used the estimation S"(a3,t) oc r~^/^ obtained from the 
following identity: 

"ti 



S°(x,ti)S^(aj,ti 



^2{t^~s)/7 
ds- ^ "^^2 



-[p{x,s)Y 



(62) 



which is obtained by substituting ti = t2 and x = x' into equation (58) using equation 
(43). 

Substituting equation (62) into equation (61), the integration of 11 is given as 
follows: 

t+At 

dt'WixA') 



d 



t+At 



dt' 



dse~^''-^y^ 



1 



7p(a;,s) 



g-2(s-si)/r 

dsi 2-fT[p{x,si)]^ 



+o(At). 

Taking the limit r ^ in equation (63) and by using the identity (48), we obtain 

"t+At 



(63) 



lim / dt'n{x, t') = -TVp{x, t)At + o(At). 
Jt 



(64) 



This evaluation shows that the coarse graining of the advection term yields the diffusion 
term in the fluctuating hydrodynamics model of Brownian dynamics. 

Finally, we obtain the change in density from time t to t + At in the limit of r ^ 
by substituting equations (57) and (64) into equation (55). The change in density to 
the order of At is written as 

dH[p] 



lim[7 {p{x, t + At)- p{x, t))] = - AtV 

T— >0 



p{x,t)V 

t+At 



6p{x,t) 

pt+i\t 

-limV-/ dt"E{x,t') + o{At). 
Jt 



TVpix,t) 



(65) 



By multiplying both sides of equation (65) by At ^ and taking the limit At 0, we 
obtain 



dp{x,t) 
dt 



7 



-p{x,t)V 



6Hv 



6p{x,t)^ 



- TVp(x, t) + ^/^Tp{x,t)C{x, t) 



(66) 



Here, we can rewrite the noise term as 
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because (60) shows that intensity in the left hand side of (67) coincides with that in the 
right hand side of (67). This is the fluctuating hydrodynamic equation for the density in 
the overdamped limit [23]. Using a technique similar to that used in the underdamped 
case, we can obtain the steady-state distribution function written as 



From the derivation of the Kawasaki-Dean formula (66), we have found that 
diffusion is caused by the advection due to a random force. In contrast, in the case 
of liquids, diffusion is caused by the liquid itself. The physical origin of the diffusion 
term obtained by our model, therefore, is different from that in a liquid system although 
the expressions of in both the cases appear similar. 

4. Discussion 

The primary objective of this study is the derivation of the underdamped nonlinear 
fluctuating hydrodynamic equation (17) along with equations (7), (13), (19), and (20). 
The starting point is the underdamped Langevin equations (1) and (2). It is a nontrivial 
fact that we obtain the closed stochastic evolution equation of the density field and 
momentum density fields from the particle description model. The exact derivation 
would have been obtained if we had taken the continuous limit using the discretization 
scheme discussed in Appendix A. The obtained evolution equation is reasonable because 
the Fokker-Planck equation obtained using our model agrees with that obtained using 
a classical liquid system, except for the form of the dissipative matrix [7]. 

The underdamped equation is unrealistic when describing the experimental 
situation of Brownian particle system. The overdamped model is more suitable as 
compared to the underdamped model, for a Brownian particle system. Nevertheless, 
the underdamped model is useful for theoretical approximations such as the mode- 
coupling theory. The underdamped model for the Brownian particles can be a basis 
for the derivation of the mode coupling equation [30]. Recently, some researchers have 
developed systematic methods for the derivation of the mode-coupling equation from the 
overdamped model for a Brownian particle system [31]. However, the derived equation 
is slightly different from the mode-coupling equation [31]. The difference might be 
eliminated if the mode-coupling equation is derived from the underdamped model. 

In addition, the underdamped model fits a liquid system. The moment of liquid 
particles should be explicitly treated in order to study the phenomena observed before 
momentum relaxation. In a liquid system, however, the field description in the nonlinear 
and underdamped cases has not been intensively studied. The present equations (21) 
and (22) for 7 = can be applied in this case. They are useful for microscopic studies 
of a liquid system. 

There are similarities and differences between our equations and the Euler equation. 
A point X in equations (21) and (22) includes not more than one particle. In 
contrast, a point x in the Euler equation includes many particles such that the 




(68) 
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thermodynamic variables are well defined. Equations (21) and (22) have been derived 
exactly, except for the condition (14). In addition, Euler equation is based on the local 
equilibrium assumption. In contrast, equations (21) and (22) can be derived without 
such assumptions. Thus, equations (21) and (22) can be used to describe the liquid that 
is not in the local equilibrium state. 

In Sec. 3.2, we have derived the Kawasaki-Dean formula by coarse-graining our 
model with respect to time under the condition (43). Similar coarse graining methods 
for the Fokker-Planck equation for the derivation of the equation of the density and 
momentum density describing a liquid system have been devised [7]. The coarse- 
graining method described in Ref. [7] does not require the condition (43). Therefore, 
the condition (43) is not required if we carry our coarse graining for the Fokker-Planck 
equation derived in Sec. 3.1. This will be investigated in our future study. 

We have found inconsistencies between the steady-state probability distributions 
(41) and (68). In Sec. 3.1, we have derived the steady-state probability distribution 
(41) for the underdamped model. We have also obtained the steady-state probability 
distribution for the overdamped model by using equation (68). The probability 
distribution (68) is not obtained by integrating equation (41) with respect to [g]. We 
guess that the inconsistencies might be related to the singularity of the delta function 
in equations (5) and (6). However, the relation between the inconsistencies and the 
singularity has not determined thus far. We will address these inconsistencies in our 
future study. 
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Appendix A. Justification of equations (14), (15) (16) and (43) 

Appendix A.l. Discretization method and justification of (14) 

In this study, we have often treated the delta function in the manner which is not 
mathematically well defined. In this section, we give the correct interpretation of these 
treatments and representations. 

First, we design a discretized cell | Aa;| that has a finite size is so small that different 
particles cannot occupy the same cell. Such a situation can be considered if the potential 
has a repulsive core within a short length r^. We denote the position of the cell by I 
introduced in Sec. 3.1. Because the cell size is sufficiently small, the map from i to J is 
an injective map. Then, the following equality is satisfied: 
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where [aj is Gauss's notation representing the maximum integer less than a and 
S[x,/Ax],ix,{t)/Ax\ = na=i,2,3^L^'i,a(i)/Ax-j,Lx,,,(t)/Ax-j- When we take the continuous hmit 
|Ax|/rj? — > 0, equation (A.l) converges to equation (14). 



Appendix A. 2. Justification of equations (15) and (16) 

The evaluation of equation (15) leads to its infinite form. This infinite form is also 
justified by discretization. First, we represent p in equation (5), in equation (6) and 
equation (10) in the discretized form as follows: 

N e 

^-^ \i\x\ 



i=l 
N 



9h 



i=l 



lAa; 



Trivially, by taking the continuous limit mentioned above, p/t and g^^ converges to 
p{x,t) and g{x,t), respectively. 

Using these discretized forms and by dividing the cell position / into two cases, we 
will prove that 

a 13 

Mf. = ^L^, (A.5) 

In the first case, consider i such that I = \_Xi{t)/Ax\. In the second case, I ^ 
[xi{t)/ Ax\ at any value of i. Equation (A.5) corresponds to equation (16) in the 
continuous limit. 

In the case of / = lxi{t)/Ax\, we can prove that the left- and right-hand sides 
of equation (A.5) are equivalent. In this case, we can calculate the left-hand side of 
equation (A.5) from equation (A. 4), so that 

, = pnt)f(t)_ (A.6) 

Since equations (A. 2) and (A. 3) reduce to pi^t = l/l^^l gf^ = pf(t)/\Ax\ 
respectively, we have 

9h9it _ [v-{t)\ ( pf{t)\ f\Ax\\ _ pmp^t) 




Equations (A. 6) and (A.7) are equivalent. 

In the case of / 7^ \_Xi{t)/Ax\ at any value of i, we prove equation (A.5) by 
defining the right-hand side by zero. From equation (A. 4) and 5/«,[zf(t)/Aa;j = 0, we 
obtain Mj^ = 0. The right-hand side of equation (A.5) is defined by zero although it 
has an infinite form because pi^t = 0. 
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Appendix A. 3. Justification of equation (43) 

The left-hand side of equation (43) is not well defined mathematically because of the 
singularity. This singularity is eliminated by the discretization of equation (43). By 
using equation (A. 2), the product of p/ j is easily calculated as 



Substituting equation (A.l) into equation (A.8) and taking the summation with respect 
to j, we obtain 



In the continuous limit, equation (A. 9) corresponds to equation (43). 
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